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ETUDE DU CHAMP ELECTROMAGNETIQUE DANS LA 
LIGNE MICRO-RUBAN BLINDEE EN UTILISANT LA 
METHODE DES DIFFERENCES FINIES 


Stefan CANTARAGIU! 


Abstract. Le présent article est le deuxiéme d'une série de trois articles publiés dans cette 
revue et qui fournit une breve présentation du livre « Micro-ondes - solution numérique », 
écrit par le méme auteur. La Méthode des Différences Finies (MDF), utilisée avec succés 
pour résoudre les problémes scalaires et vectoriels les plus difficiles de l'électrodynamique, 
facilite l'approximation des équations de Helmholtz dans un nombre fini du domaine 
analysé. Les avantages de la méthode résident dans la possibilité de réunir les deux 
domaines en un seul probleme et dans sa mise en ceuvre plus rapide, ce qui permet de 
définir facilement les conditions aux frontiéres. Dans la derniére section de l'article, les 
résultats et les graphiques des composantes du champ électromagnétique sont présentés. 
Les inférences de simulation a la frontiére des domaines air et diélectrique sont également 
mises en évidence. Les modes de propagation hybrides d'ordre supérieur sont déterminés 
en Sélectionnant les valeurs propres une a une par ordre croissant de fréquence. 


Keywords: Méthode des Différences Finies (MDF), Méthode des éléments finis (MEF), 
équations de Helmholtz, valeurs propres et vecteurs propres, conditions aux frontiéres, 
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I. Introduction 


L'article vise l'étude du champ électromagnétique dans la ligne micro-ruban 
blindée en utilisant MDF, qui est l'une des méthodes numériques les plus 
performantes, qui a été utilisée avec succés pour résoudre les problémes les plus 
complexes de l'électrodynamique, et qui refléte le comportement du champ 
électromagnétique dans le structure analysée. 


Compte tenu des techniques d'approximation numérique, utilisées pour résoudre 
des équations aux dérivées partielles, les travaux publiés dans la littérature 
spécialisée [1]+[10], qui traitent de l'analyse des modes de propagation hybrides et 
des propriétés de dispersion des lignes micro-ruban, peuvent étre divisés en deux 
groupes. Le premier groupe concerne des méthodes numériques précédées d'un 
traitement analytique important, tandis que le deuxiéme groupe se caractérise 
plutét par un traitement analytique extrémement rudimentaire, toute la difficulté 
étant transférée aux procédures de calcul disponibles sur le marché. 


' Chercheur principal, membre correspondant de I'Académie des Scientifiques Roumains. 


16 Stefan CANTARAGIU 


Les méthodes numériques font principalement référence a la FDM et a la méthode 
des éléments finis (MEF) pour résoudre les problemes électrodynamiques les plus 
complexes, tels que le comportement du champ électromagnétique dans la ligne 
micro-ruban blindée et mettent en évidence les frontiéres des formules utilisées 
par la méthode quasi-statique aux fréquences micro-ondes plus élevées. 


MDF offre des avantages significatifs en cas d'analyse de_ structures 
électrodynamiques; ces avantages seront mis en évidence dans le contenu de 
l'article. 


II. Formulation 


MDF est appliqué pour approximer les solutions des doubles équations de 
Helmholtz: 


2 c 2 : 
O*H,(xy) , O*Hz(xy) +k5H,(x,y)=0, 


Ox? ay? 
a7EZ(x%,y) , O7EZ (xy) , 22 = 
ax2 + ay2 +ksE, (x, y) =f), 


dans le cas des problemes de Neuman et Dirichlet (k représente le nombre d'onde 
longitudinal, et l'indice 6, avec d=1+2, met en évidence les domaines air et 
diélectrique). 


Les solutions des équations de Helmholtz sont déterminées, continues sur le 
domaine compact D UT (I est la fronti¢re du domaine D), délimité par la section 
transversale par la ligne micro-ruban blindée (figure 1) et qui satisfont aux 
conditions aux frontiéres: 


oH = R 
ry = 0, nest le verseur normal a I. 
Fr 


E.\|r — 0. 


Pour réduire le probleme 4 un probleme numérique, le maillage de la figure | est 
considéré: 
X = Xo + JAx, 


y=yo + iAy, ijEZ, 


1 ‘2 
qui couvre le domaine pG) Ur () (pour des raisons de symétrie, seule la moitié 
de la section de la ligne micro-ruban est utilisée). 
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‘a 1 
Les points de pG) Ur @), qui appartiennent au réseau et qui ont au moins un 


0) 


point voisin extérieur au domaine forment la frontiére du réseau, 


Y‘p | Ax 


Q 


Figure 1. Le réseau des cellules rectangulaires des cellules associées a la cellule 
élémentaire 


Dans chaque noeud (j,i) (qui correspond au point de coordonnées 
Zz 


(xq + jAx, Yo + iAy) = (x;,y;)) dans le domaine D,*’, la premiére équation de 
Helmholtz est approximée par une équation aux différences finies qui utilisera les 
approximations des dérivées de la composante longitudinale du champ 
magnétique (la notation @ = H, est utilisée): 


9g\  , $jt1 iP j-i 
aR = 2Ax : (la) 
O7o\ _ Dj41i-2Ojit Ojai 
(SS di fia Ax2 ’ (1b) 
OG\ _, Pjiti—~Pj ia 
Go) ~ 2Ay : (Ic) 


(<2) ™ Pjits 20 jit Oji-a (1d) 


dy? ji Ay 


Les dérivées de premiére ordre (1a) et (1c) sont utilisées afin d’établir les relations 
entre les composantes longitudinales du champ magnétique, a la proximité de 


Vécran électrique. La notation 7 = “a Es est utilisée afin de simplifier les 


relations entre les composantes tangentielles 4 la surface de séparation des deux 


18 Stefan CANTARAGIU 


milieux. Les approximations par différences finies des dérivées de la composante 
longitudinale du champ électrique sont déterminées de la maniére suivante: 


L(o) = 4° + k§p = 0, (2a) 
L@) = 4*n + kZn = 0. (2b) 
Utilisant les notations: 
Ax k2 
A= kAx?, R == et T= 2 (3) 


les équations différentielles deviennent: 
Ag; = 2(1 + R*)0;,i — 641,35 — Oj-1i — R7 Gj 541 — R° Oi, 4 
Anji = 2(1 + R*)nji — Nj — jai — R?Nj isa — R?Nj ia. (5) 
pour I’air et: 
Ag; ; = 2t(1 + R7) Gj; — TH 413 — THj-15 — TR? Gj 541 —TR*Gji-4, (© 
Anj¢ = 2t(1 + R?)nji — Tj4a¢ — THj-1g — TR? Njina —TR?Nji-1, (7) 
pour le milieu diélectrique. 


L’introduction du paramétre 7 est trés importante, car elle nous permet de 
déterminer la constante de temps du champ électromagnétique dans la ligne 
micro-ruban (les deux domaines sont réunis dans un seul probléme). 
5 ar ae = 
Pour chaque point de l’ensemble D,~ une équation aux différences finies est 
écrite, similaire aux (2.7) et (2.9), pour la composante @ et aux (2.8) et (2.10) pour 
la composante 7. Parce que les fonctions inconnues @ et 7 sont données sur la 
1 
courbe /} les procédures de transfert des conditions a la frontiére sur r@) (qui 
n’est plus une courbe mais un ensemble de points) sont utilisées afin d’obtenir les 


r () 


dépend des valeurs des points voisins de N dans domaine D 


conditions sur l’ensemble 
ih 
point NE r@) 
sur I’. 


. La valeur attribuée aux fonctions @ et 7 dans un 


() 


et 
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III. Equations aux différences finies sur les frontiéres 


A la proximité d’écran électrique, oti les conditions a la fronti¢re imposent [8]: 
ag 
; ;=0 et (55) =0 
jai on je 


et, a la proximité d’écran magnétique (situé dans le plan x =0 sur la figure 1), ot 
les conditions a la frontiére sont [8]: 


—— a) - 
ji = 0 et ( Pa 


les équations (4)+(7) se modifient. Sur la fronti¢re PQ (voir la figure 1), ot: 


OO; 
oy 


= Oetn; ;=0, 


les équations (1c) et (2d), au-dessus de la fronti¢re PQ, permettent écrire pour le 


@), 


noeud fictif (j,i—1), placé a l’extérieur du domaine D? 


. _ 4 fictif 
Dj i+ ~_ jini 

_ _, fictif 
1j,i+1 — “Nj i-41 


La derniére équation concernant l’asymétrie relative a l’écran électrique de la 
A A Z 5 QWEg agi a 
composante longitudinale du champ électrique (7 = “p Fz) sera utilisée afin de 


déterminer les composantes transversales du champ électrique. Compte tenant de 
la premiere équation ci-dessus, |’équation (4) devient: 


Ad; ; = 2(1 + R*) 0; i — 6415 — Ojai — 2R7Gj 544 (8) 


De maniére similaire, sur la frontiére SR, le noeud fictif (j,i+1), placé a 
1 


po8 : 2 ne sys 2 
Vextérieur du domaine D, , sous la frontiére SR, est utilisé et on trouve: 


fictif _ 4 | 
j.it1 = $j i-1, 


fictif_ 
Nita — —1ji-1 
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Dans ces conditions, |’équation (6) devient: 


Aj; = 2t(1 + R7)O;; — Thj413 — THj-1,; — 2TR?G; 3-4 (9) 
adj; 
Sur la fronti¢re QR (figure 1), ot “Pit — 0 et ji = 0, compte tenant des 


équations (la) et (2b), il résulte: 


fictif _ .~ | 
j+1,i = $j-1 


fictif 
Tpeag ~~ 8j-1F- 


1 
Dans ce cas, le nceud fictif est placé a l’extérieur du domaine DS , dans la partie 
droite de la frontiére. Par conséquence, |’équation (4) devient: 


Ad; ; = 2(1 + R*) ji — 263-15 — R76; 541 — R? Oj 4-1, (10) 
et l’équation (6) devient: 


Aj; = 2t(1 + R?)G;; — 27,4; — TRG; 4 — TR? Gj 5-1 (11) 


On; ae 
Sur la frontiére PS (figure 1), o1 #* = 0 et $;,; = 0, il résulte: 


Ox 
_ _ fictif 
1j+1,i = Vj-1i 
= fictif 
Pj +1 — P-L 


Le noeud (j—1,i) est lui aussi un noeud fictif placé en dehors de la cellule 
élémentaire, 4 gauche de la fronti¢re PS. Dans ce cas, |’équation (5) devient: 


Anj¢ = 2(1 + R*)nji — 2nja1i — R?Nj i441 — R?nji-1, (12) 

et l’équation (7) devient: 
An jz = 2t(1 + R?)nji — 20 j41,; — TR? j¢44 — TR? 5-1. (13) 
L’antisymétrie de la composante longitudinale du champ magnétique (@ = H-) 


relative a l’écran électrique est utilis¢e afin de déterminer les composantes 
transversales du champ magnétique. A la fronti¢re de séparation des domaines, 
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définies par le produit cartésien {y = y,} X FE. ‘|. compte tenant de la loi de 
Gauss pour le champ magnétique (1d), la condition suivante est accomplie: 


7B, —7iBp —0 (14) 


Les indices utilisés dans |’équation (14) mettent en évidence la densité du flux 
magnétique dans |’air (le premier) et respectivement la densité du flux magnétique 
dans le milieu diélectrique de la sous-couche (le deuxieme). Parce que 
Ht = Hz = 1, on peut écrire: 


wa 


Hy =Hyp|_» x € (5,5) (15) 


De plus, selon la relation entre les composantes tangentielles des champs 
électrique et magnétique a la frontiére de séparation entre les domaines air et 
diélectrique, respectivement H,9, exprimé sous forme de séries rapidement 
convergentes utilisant des systemes orthogonaux de polynémes et des fonctions de 
Chebyshev, il en résulte: 


wa 
Considérant les relations entre les composantes tangentielles du champ 


électromagnétique a la surface de séparation entre les domaines air et diélectrique, 
les relations (15) et (16) deviennent: 


5G, , 91, _ (_ 96a anp 

— Fa 4 a= 7 ( ay + Er2 i) (17) 
O¢4 , ON _ 9¢p Onp 
ax + Oy 7 r( Ox t Er2 me) (18) 


Afin d’établir la liaison entre les composantes des champs électrique et 
magnétique, tangentielles a la surface de séparation des milieux, la notation P41 
est introduite; elle met en valeur les valeurs des composantes longitudinales du 
champ magnétique, correspondantes aux points (j,i-/), placés dans l’air; aussi, la 
notation $; i+; représente les valeurs correspondantes aux points (j,i+1), placé 
dans le milieu diélectrique. Les équations aux différences finies correspondantes 
aux points situés sur la surface de séparation des milieux peuvent étre écrites une 
fois utilisant la relation (4), quand I’interface est abordée du cété air: 


Adj; = 2(1 + R7) 9; i — Oj41i — O)-15 — R? Oj i401 —P?OA4, (19) 
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et la deuxiéme fois en utilisant la relation (6), quand l’interface est abordée du 
cété diélectrique: 


Ag; ; = 2t(1 + R*)Oj 5 — THj415 — THj-1,i — TR°DP ig — TR?Hji-1 (20) 


Afin d’éliminer les valeurs inconnues PFi-4 et ?; j+1 dans les équations (19) et 


(20), les conditions (17) et (18) sont écrites sous une forme utilisant les 
différences finies: 


D A 

Pj+1i—Pj-1i Npita ja Njita Nj ia 

Tt — 1 a =. Ee. Tt et 
( ) 2Ax r2 2Ay 2Ay (21) 

A D 
Aj+ai—Nj—-1i Pj ita Pj i-1 — Pjita~Pji-n 
E->T —1 = E,9T 

( T2 ) 2Ax r2 2Ay 2Ay : (22) 


Suivant |’élimination de Pi et de PPivt dans |’équation 22), et utilisant les 
relations (19)+(20), il résulte: 


1 
Agji = (1 + T)(1 + R7)Gj5 — 5 (1 + 1) Oj425 — R° Oj i+ — 
1 1 1 
—5 (1 + 1) bj-1,5 — TRG; 5-1 — FRCL — &p2t) jars + [RO — &r2t)nj-1i- (23) 


Les valeurs des composantes longitudinales du champ magnétique sur la ligne 
métallique (notée par VU sur la figure 1) sont déterminées en utilisant la relation 
(8) pour le cas ot le réseau de droites aborde la ligne du cété air et en utilisant 
l’équation (9), quand celle-ci est abordée du cété diélectrique. 


L’équation correspondante aux composantes longitudinales du champ électrique 
est obtenue de maniére similaire en utilisant la relation (21), partant des équations 
aux différences finies correspondantes aux points situés sur la surface de 
séparation des milieux, écrites par l’interméde de |’ équation (5) (quand |’ interface 
est abordée du cété air) et par l’interméde de |’équation (7) (quand I’interface est 
abordée du cété diélectrique): 


1+ ,2T 1+ €,2T 2R? 
7 1=2(7-)a R? u-G*) _ 
1) jit 14 Er? ( + )n ji ie E+2 I jt1t bs 41 


*r2t) 2e,2TR? ( t-1 ) ( 1 ) 
= He ye Sep en ee a 
cen Ij-1i Er2t1 1ji-1 1+e,2 Pj+1i 1te,2 $j-1i (24) 


1 j,i+1 


Les équations (23) et (24) sont valables au plus grand nombre de points sur la 
surface de séparation entre les domaines air et diélectrique. 
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Il existe deux points sur la surface de séparation entre les deux environnements 
présentant des particularités, qui seront analysés plus en détail. 


L'un d'eux est le point T de la figure 1, ou l'interface rencontre la frontiére QR, et 
la relation (24) ne peut pas étre utilisée, car ici, n i= 0, respectivement: 


ani 
Np = By? = 0, 


et en substituant l'identité ci-dessus dans l'équation de Helmholtz (2b), il en 
résulte: 


077; 5 
—=0 22 
La deuxiéme condition au voisinage du mur électrique impose que: 


995i 


Ox = 0. 


Par conséquent, on obtient, conformément aux relations (1a), (2d) et (2b): 


Pja1i = 9-13. 
pitt — “Nj i-1> 
j+1i — —Nj-1j 


Ces conditions déterminent le changement de la relation (23), comme suit: 


Adj; = (1 + 7)(1 + R*)G; is — 1 + Gj 41,5 — R? Oj 541 — 
—tR?p; «1 — R(1 — &,27) 41, (26) 


L'analyse du comportement des composantes longitudinales des champs électrique 
et magnétique au voisinage de l'autre point de la surface séparant les domaines, 
noté U sur la figure 1, se fait de maniére analogue. 


IV. Probléme de valeurs propres 


Si les Equations @; ; et ;; sont sélectionnées a chaque fois parmi les relations 
(4)+(26), qui particularisent chaque noeud séparément du maillage représenté a la 
figure 1, on obtient un probleme de valeurs propres et de vecteurs propres: 


AX =A,X, (27) 
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ou 


X est un vecteur propre; 

A; (k représente le nombre d’équations aux différences finies dans le systéme) 
sont les valeurs propres; 

A est une matrice carrée, dont |’étude sera reprise quand les logiciels de calcul 
seront détaillés. 

Le vecteur X est formé par la réunion des fonctions @ ji et Nj i, C’est-a-dire: 


xy P11 

X2 P12 

X3 : 
: Pim 
: a1 
. = Darn ! k=nxX m+sx p 
. 11 
: 112 
. Nip 
: 21 
k 


Isp 


La réunion des composantes longitudinales des champs magnétique et électrique 
dans un seul vecteur est imposée par le fait que les équations de continuité a la 
surface de séparation des milieux, (23) et (24), utilisent les composantes 
longitudinales du champ électrique mais aussi celles du champ magnétique; cette 
réunion est possible grace au paramétre T . 


La solution du probléme matricielle (27) permet a déterminer les composantes 
longitudinales puis transversale des champs électrique et magnétique. 


V. Simulation informatique 


Il s'agit d'une véritable structure de ligne micro-ruban blindée, fréquemment 
utilisée dans le domaine des radars, travaillant dans la gamme des ondes 
centimétriques et millimétriques, pour laquelle w = 1 mm, a = 3,5 mm, 
y, = 0,5 mm, yz = 2mmet £,7 = 9. 
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L'analyse du comportement du champ électromagnétique au voisinage des parois 
électriques et magnétiques est plus facile, si les nceuds du maillage qui recouvre la 


cellule élémentaire sont numérotés (i=/~+/0, j=1<8, Ax = =0,25 mm, et Ay = 


14 
"2=0,25 mm). 


Le maillage se compose de huit colonnes et dix lignes. Lorsqu’un maillage plus 
fin est utilisé, les erreurs introduites par les approximations FDM sont plus petites. 


a eee fas 3 a 
Les noeuds situés a I'interface entre les deux milieux, ot x € [o, | (on note que la 


zone ou se situe la bande est également analysée), sont numérotés deux fois, car 
ici deux lignes (i=7 et 1=8) sont utilisées, se chevauchant pratiquement. 


Méthodes de détermination des composantes longitudinales du champ 
électromagnétique et des parameétres de propagation 


La détermination des composantes longitudinales du champ électromagnétique 


E.(n) et H-(@) est réalisé a l'aide de l'environnement de développement intégré 
Matlab. 


Le tableau 1 ci-dessous résume les conditions aux fronti¢res auxquelles sont 
soumises les composantes longitudinales des champs électriques et magnétiques, 
conformément aux noeuds de la figure 1. Ces conditions sont prises en compte lors 
de l'élaboration du programme de calcul. 


Le probléme aux valeurs propres, généré apres écriture des équations aux 
différences finies, correspondant aux composantes longitudinales du champ 
électrique et magnétique, calculées pour chaque nceud du maillage, s'écrit sous la 
forme: 


(Ad — A)X = 0, 


ou I est la matrice unitaire, A est la matrice carrée définie par (27) ; dans le cas de 
l'exemple du chapitre, la taille de la matrice comporte 12.100 éléments — 110 x 
110; k représente le nombre d'équations aux différences finies (k=1+110). Pour 
réduire l'ordre de grandeur des erreurs qui apparaissent en remplacant les 
équations de Helmholtz par des équations aux différences finies, le nombre 
d'éléments pourrait étre augmenté en conséquence. 


La numérotation des noeuds du maillage, qui correspondent aux valeurs des 
composantes longitudinales du champ, peut étre choisie, en principe, par choix - 
dans le cas présent avec les noeuds qui correspondent aux valeurs des composantes 
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longitudinales du champ, champ électrique, situé dans la partie supérieure gauche 
de la cellule élémentaire. Le numéro de séquence du neeud est placé dans la partie 
inférieure droite de celui-ci. 


Tableau 1 


PQ et SR nOnZOuUE 
vertical 
vertical 


en dessous de la bande VU horizontaux 
point U vertical 


les noeuds sur la surface de air - diélectrique | nf; = n> 07; = 07; 
séparation entre les supports Hya = Hyp; Hyg = Hyp 
diélectriques 


Apres avoir épuisé les valeurs des composantes longitudinales, les noeuds sont 
renumérotés, cette fois pour mettre en évidence les valeurs des composantes 
longitudinales du champ magnétique aux points du maillage, en partant de la 
partie inférieure gauche de la cellule élémentaire. Le numéro de séquence du 
noeud est placé dans la partie supérieure droite de celui-ci. 


Chaque nceud du maillage a été numéroté deux fois, la premicre fois pour mettre 
en évidence la valeur de la composante longitudinale du champ électrique, et la 
deuxiéme fois pour mettre en évidence la valeur de la composante longitudinale 
du champ magnétique. Les noeuds ot les valeurs des composantes longitudinales 
des champs électriques et magnétiques sont nulles recoivent le numéro 0. 


Pour déterminer les valeurs des composantes longitudinales qui correspondent aux 
noeuds sur les frontiéres, selon les étapes présentées, des séries de noeuds factices 


sont utilisées. Ces derniers sont situés en dehors du domaine p& 2) , et la saisie des 
valeurs qui leur correspondent dans la matrice A se fait en thant compte des 
relations (8)+(13). L'antisymétrie des composantes longitudinales par rapport aux 
frontiéres est indiquée sur le dessin avec l'appui du signe négatif, et leur symétrie 
du signe positif. 
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A la surface de séparation entre les deux domaines, définie par la relation 
cartésienne {y=y1} x[E, =) les relations (23)+(26) sont utilisées. Le probleme 
des valeurs propres et des vecteurs propres est généré en écrivant dans leur ordre 
naturel les 110 équations correspondant a chaque noeud du maillage (i=1+10, 
j=1+8), de 1 a 110, en sélectionnant chacune d'entre elles, parmi les relations 
(4)+(26), l’équation qui singularise chaque noeud. L'introduction de la matrice A 
dans le programme de calcul devient plus facile par partitionnement. On considére 
l'observation que les éléments non nuls de la matrice sont regroupés dans sa zone 
diagonale (fait da a l'utilisation de la procédure "d'approximation croisée" lors de 
l'application du MDF). 


La procédure de calcul est la suivante : pour la valeur z choisie par le programme 
parmi les 110 valeurs propres, on détermine la valeur propre correspondant a la 
fréquence de travail la plus basse de la gamme (pour t= - 5 la valeur propre, qui 
désigne le mode hybride fondamental de propagation depuis la ligne, correspond a 
la constante de phase f=1131,62 rad/m et a la fréquence f7=19,509 GHz). 


La valeur propre avec la plus petite valeur négative désigne le mode de 
propagation fondamental hybride dans la ligne de transmission micro-ruban 
blindée. 


a Elvin} a of [Him] 


Figure 2. La variation des composantes et correspondant au mode hybride fondamental 
de propagation dans les neuds du maillage ((x;,yi), j=1+8 et i=1~+10). 


L'étape suivante est la détermination du vecteur propre a partir de la structure 
duquel les composantes longitudinales du champ électrique sont séparées de celles 
du champ magnétique. 


La figure 2 montre les graphiques des composantes électriques et magnétiques 
longitudinales qui mettent en évidence le fait que les variations majeures se 
situent, comme prévu, dans la zone d'interface entre les deux domaines et dans la 
zone diélectrique du substrat. Les modes de propagation hybrides d'ordre 
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supérieur sont déterminés en sélectionnant les valeurs propres une par une par 
ordre croissant de fréquence. 


Les figures 3 et 4 montrent les composantes longitudinales du champ 
électromagnétique, déterminées avec les deux valeurs suivantes des constantes de 
phase, choisies par ordre croissant de fréquence (f2=60,385 GHz et fs=148,9807 
GHz)) et qui correspondent a la supérieure modes de propagation hybrides d'ordre 
2. et 3. 


405 E23 [V/m] 


4 
x [mm] 


Figure 3. Variation des composantes longitudinales du champ électrique Ezz et Ez3 et 
dans les neeuds du maillage ((x;,yi), j=1+8, et i=1~+10) 


085 H, 3[A/m] 


Figure 4. Variation des composantes longitudinales du champ magnétique Hz et Hz3 
dans les neeuds du maillage ((x;yi), j=1+8, et i=1~+10) 


Méthodes de détermination des composantes transversales du champ 
électromagnétique 


Les composantes transversales des champs électrique et magnétique (E,, Ey, H, 
et Hy) sont calculées dans lequel interviennent les dérivées des composantes 


1 1 
longitudinales. Pour un point du domaine pG) Ur () le calcul des dérivées se 
fait a l'aide des relations (1a), (1c), (2a) et (2c). Les dérivées correspondant aux 
points sur les frontieres (parois électriques et magnétiques) sont calculées en 


1 
introduisant quelques points fictifs, placés en dehors du domaine p&@), 
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Dans le cas du mur électrique vertical, délimité par la frontiére QR, ol —— wees = Oet 


1g i=0, en considérant les relations (1a) et (2b), il s'ensuit: 


= =o, - (28) 
bit = —n7;, i=1+10 (29) 
954 


Pour les murs électriques horizontaux, ot i. = 0 et ;;=0, en considérant les 


relations (1c) et (2d), il s'ensuit : 


- correspondant a la frontiére PQ 
%j2=)5 7 (30) 
nj2 = npg! j=128; (31) 


- correspondant a la frontiére SR 


on = d9 (32) 

ni it — ng, j=128 (33) 

Pour le mur magnétique, noté PS, ot aes = 0 et @ji=0, en considérant les 
relations (2a) et (1b), il en résulte: 

Mi = Mai (34) 

bo; = — 05°", i=1=10 (35) 


Pour le calcul des dérivées des composantes longitudinales a la surface de 
séparation entre les deux domaines, définie par la relation cartésienne 


ty=yu)x[F, | (frontiére UT), les noeuds fictifs (j, 6) et (j, 9) sont introduits, que 


l'on retrouve dans les relations correspondant au champ magnétique: 


i-7> (24) — tisFjs (36) 
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i=8=>( (37) 


as) _ Pjs~9j7 
ay j8 2Ay ° 


ou j=428. 


L'élimination des valeurs correspondant aux noeuds fictifs se fait a l'aide de la 
relation (4), écrite pour i=7: 


Adj 7 = 2(1 + R?)b;7 — bj417 — 0-17 — R° jg — R° O76 (38) 
et de la relation (6), écrite pour i=8: 
Agjg = 2t(1 + R*)bjg — thj+1g —THj-18 — TR’ Pg — TRG; 7 (39) 


De maniére analogue, les dérivées des composantes longitudinales du champ 
électrique sont calculées, respectivement pour i=7 et pour i=8 pour j=4+8. 


L'élimination des valeurs correspondant aux noeuds fictifs se fait, cette fois a l'aide 
de la relation (5), écrite pour i=7: 


Anj7 = 2(1 + R?)nj7 — nj+17 — Nj-1,8 — R?njg — R?nZe (40) 
et de la relation (7), écrite pour i=8: 
Anjg = 2t(1 + R*)njg — TNj+18 —™Mj-1g—TR*nPy—TR*Nj7 (41) 
De maniére analogue, on calcule les dérivées des composantes longitudinales des 


champs électriques et magnétiques, correspondant aux noeuds situés sur la bande 
VU, ow: 


fictif = 
i<7=> oe a (42) 
fictif __ 
Nig =—Nje 
et 
plies =>; 
acs lag (43) 
fictif _ 
n>) =—Njs 
pour j=/+3. 


L'élimination des valeurs des composantes du champ magnétique, correspondant 
aux noeuds fictifs, se fait a l'aide de la relation (4), écrite pour 1=7 et de la relation 
(6), écrite pour i=8. 
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Les figures 5 et 6 montrent les graphiques des composantes transverses des 
champs électrique et magnétique, correspondant au mode hybride fondamental de 
propagation en chaque noeud du maillage. 


Conformément a la méthode d'identification des modes de propagation d'ordre 
supérieur, décrite dans la section IV, les figures 7 et 8 montrent les composantes 
transversales du champ électromagnétique, déterminées avec les deux valeurs 
suivantes des constantes de phase, choisies par ordre croissant de fréquence et qui 
correspondent a des modes hybrides supérieurs d’ordre 2 et 3 pour Ey et Ey 3, 
E,.2 et Ex 3. 


Vin) 
f ops 


Figure 5. Variation des composantes E), et Ey correspondant au mode de propagation 
hybride fondamental dans les neeuds du maillage (x;,yi), j=1+8 and i=1+10 


fy, x(n) 


Figure 6. Variation des composantes Hy, et Ay correspondant au mode de propagation 
hybride fondamental dans les neeuds du maillage (x;,yi), j=1+8 and i=1~+10 
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200 Ey 3 [V/m] 
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x [mm] 


Figure 7. Variation des composantes transverses du champ électrique Ey. et Ex.3 dans 
les neuds du maillage ((x;,yi), j=1=8, et i=1~+10) 


200) Ey.2 [Vim] 1805 Ey 3 [Vén] 


Figure 8. Variation des composantes transverses du champ électrique Ey et Ey 3 dans 
les neuds du maillage ((xj,yi), j=1=+8, et i=1~+10) 


Variation des composantes transversales du champ électrique Hy 2, Hy3, Hy» et 
Hy 3 sont calculées de la méme maniére. 


VI. Remarques finales 


MDF est une technique numérique particuli¢rement puissante, utilisée avec succés 
dans la résolution des problemes scalaires et vectoriels les plus difficiles de 
l'électrodynamique [8]+[10] et [14], qui permet l'approximation des équations de 
Helmholtz avec des différences finies en un nombre fini de points de le domaine 
analysé. 


L'un des avantages du MDF réside dans la possibilité de réunir les deux domaines 
analysés dans un seul probléme et, implicitement, dans la mise en ceuvre et la 
résolution plus rapides du modeéle mathématique. Dans le méme temps, le MDF 
facilite facilement les problémes de manipulation des frontiéres dans le but de 
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définir les conditions aux frontiéres, respectivement au voisinage de la bande ou a 
l'interférence entre les domaines air et diélectrique. 


A cet égard, l'étude choisit de satisfaire les conditions imposées au champ 
électromagnétique a l'interférence entre les domaines air et diélectrique, et il est 
ainsi possible de rassembler les composantes longitudinales des champs électrique 
et magnétique au sein d'un méme probléme de valeurs propres et vecteurs. Suite a 
l'identification des vecteurs propres, les valeurs des composantes longitudinales 
puis transversales des champs électriques et magnétiques sont déterminées en tous 
points du maillage dans les domaines analysés. 


L'ordre de grandeur des erreurs qui apparaissent en remplacant les équations de 
Helmholtz par des équations aux différences finies est égal a l'ordre de grandeur 
du carré du pas de maillage sélectionné ((Ax)* et (Ay)7), qui doit étre le plus petit 
et le plus fin possible. 
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